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ABSTRACT 

We report on the initial analysis of HerschelfHIFl carbon monoxide (CO) observations of the Martian atmosphere performed between 1 1 and 16 
April 2010. We selected the (7-6) rotational transitions of the isotopes 13 CO at 771 GHz and C I8 and 768 GHz in order to retrieve the mean 
vertical profile of temperature and the mean volume mixing ratio of carbon monoxide. The derived temperature profile agrees within less than 
5 K with general circulation model (GCM) predictions up to an altitude of 45 km, however, show about 12-15 K lower values at 60 km. The CO 
mixing ratio was determined as 980 ±150 ppm, in agreement with the 900 ppm derived from ifer.sc/ie//SPIRE observations in November 2009. 
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1. Introduction 

Carbon monoxide was first definitely detected spectroscopically 
in the Martian atmosphere by Kaplan et al. (1969). From obser- 
vations performed in 1967 between Mars L s - 110° and 133° a 
volume mixing ratio (vmr) of 800 ppm was deduced Kakar et al. 
( 1977 ) reported the first microwave detection of CO performed 
in 1975, when observing the J = 1-0 rotational transition during 
L s - 339°. They found a vmr of 1900 ppm assuming a disk- 
averaged continuum brightness temperature of 200 K and a con- 
stant height profile. Good & Schloerb ( 1981 ) observed the same 
transition in 1980 during L s = 86° and found 3200 ppm. The re- 
analysis of these data by Clancy & Muhleman ( 1983) resulted in 
quite different values: 800±400 ppm from the Ka kar etal.|( 1977|l 



observation and 1400 ± 500 ppm from the |Good & Schloerb 



(JT98Tb observation by assuming a surface brightness tempera- 
ture of 206 K, a ±10 K variation in model atmosphere tempera- 
tures, and a surface pressure of 6.5 hPa. They found 2200 ± 1000 
ppm, based on their own observations of the CO (2-1) line per- 
formed in 1982 (L s = 75°) with the same model assumptions. 
Clan cy et al.| (1990) derived a mixing ratio of 600 ± 150 ppm 
from 1988/9 observations (L s = 197° - 341°) of the 12 CO and 
13 CO transitions (1-0 and 2-1) and surface brightness tempera- 
tures between 192 and 216 K. They did not find any significant 
temporal variability of CO. [Rosenqvist et al. ( 1990) found a CO 
vmrs between 800 and 1200 ppm depending on the orography 
from Phobos/ISM observations. Lellouch et al. (1991) provided 
for the first time a mm-wave map of CO (1988, L s = 270°) and 



deduced no significant variation from their derived value of 800 
± 200 ppm over the Mars disk. While the 230 GHz observa- 
tions provided reasonable surface brightness temperatures, the 
values in the 1 15 GHz band are too high (270 K), very likely 
due to sideband ratio calibration errors. Billeb aud et al.| (l998) 
retrieved values between 450 and 1 150 ppm from ground-based 
IR observations in 1990/1. Also from ground-based IR obser- 
vations, Krasnopolsky (2003) found a hemispheric asymmetry 
with the CO vmr increasing from 850 ppm at 23° N to 1250 ppm 
at 50° S (at L s - 1 12°). Follow-up observations of Krasnopolsky 
(2007) between 80° N to 75° S over 4 seasons show even larger 
variations between 700 and 1600 ppm. Encrenaz et al. (2006) re- 
port about a seasonal variation in CO over Hellas by a factor of 2 
from 2.3 //m CO band observations by OMEGA/MEX. Infrared 
observations of PFS/MEX result in similar variabilities over lat- 



* Herschel is an ESA space observatory with science instruments 
provided by European-led Principal Investigator consortia and with im- 
portant participation from NASA. 



itude and season, i.e. with a mean value of 1 100 ppm (Billebaud 
|et al.|2009) . Finally, [Smith et al.| ( p009) ) derived seasonally and 
globally averaged 700 ppm CO vmr from CRISM/MRO obser- 
vations, but with strong seasonal variations at high latitudes. The 
summertime near-polar CO mixing ratio was observed to fall to 
200 ppm in the south and 400 ppm in the north during the time 
the carbon dioxide was sublimating from the seasonal ice caps. 

The determination of the vertical profile of temperature and 
CO (together with studies of the carbon and oxygen isotopic 
ratios in CO) over different solar longitudes of Mars is one of 
the goals of the Herschel key programme "Water and related 
chemistry in the solar system" (Harto gh et al.||20"0"9) >. Here we 
present first results of observations performed in April 2010 
around L s = 76°. In this paper we derive the CO mixing ra- 
tio and temperature profile from HIFI observations. The latter is 
also required for analysing other species not treated in this paper. 
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For simplification we analyse two CO lines (the 7-6 transitions 
of the isotopes 13 CO and C I8 0) observed at the same time and 
in the same sideband. 



2. Herschel/H\F\ observations 

All observations were carried out in HIFI's dual beam switch 
mode (de Graauw et al. 2010 Roelfsema et al. 2010 1 on op- 



erational days (OD) 332, 333, 334, and 337, corresponding to 
1 1-16 April 2010 or L s = 75.8° to 78°. Dedicated line observa- 
tions include CO (5-4, 6-5, 8-7), 13 CO (7-6), C n O (7-6), and 
C ls O (7-6). HIFI spectrometers can resolve line shapes with a 
spectral resolution of 140 KHz (High Resolution Spectrometer; 
HRS) or 1.1 MHz (Wide Band Spectrometer; WBS). The inte- 
gration times range from 93 s for the strongest I2 CO lines to 
9289 s for the C n O observation. Numerous other CO-lines were 
detected in the band scans or as side products of other dedi- 
cated line observations, appearing in the intermediate frequency 
either in the same or other sideband (double sideband conver- 
sion), some of them with excellent signal-to-noise ratios. The 
first set of data was available about a week after the observations 
and was processed with the standard HIPE v3.0.1 modules ( |Ott| 
2010) up to level 2. This data set was not complete yet; for in- 
stance, the data of the HRS was only partly available, and point- 
ing products therein had no entries. The flux calibration by the 
Herschel Science Centre (HSC) was still in progress and some 
calibration errors occurred in this first data set. Furthermore, it 
turned out in the band scan data that the line amplitudes of the 
detected CO lines were not always exactly the same in the upper 
and lower sidebands. Last but not least, the HIFI spectra suf- 
fered from a relatively large baseline ripple (Fig.fTb. Fortunately, 
the frequency of the ripple generally is the same in the vertical 
and horizontal polarizations; however, phase and amplitudes are 
different so they need a dedicated treatment before both polar- 
izations can be averaged. Solutions of these problems are un- 
der way, but only partly available now. Therefore we decided 
to analyse the data with the following boundary conditions: (i) 
line-to-continuum ratio rather than absolute fluxes (requires as- 
sumptions on the surface brightness), (ii) only horizontal polar- 
ization (hence decreases the signal-to-noise ratio), (iii) only lines 
in the same sideband, (iv) only WBS spectra. The surface pres- 
sure averaged over the visible Mars disk changes diurnally in 
the order of 10% (derived from the European Mars Climate Data 
Base, EMCD v4.1,|Forget et al.|1999t|Lewis et al.|1999|and our 



own general circulation model, IHartogh et al. 2005 ; Medvedev 
& Hartogh 2007). Therefore we introduced boundary condition 
(v) to analyse only lines observed at the same time. The only 
data that fulfil these criteria are the (7-6) transitions of 13 CO 
and C ls O, which were observed on 13 April 2010 at 05:33 UT 
(Obs. Id. 1342194686). The integration time for this observation 
was 900 s. 



3. Analysis and discussion 



Compared to cometary observations of HIFI (Hartogh et al. 
2010, de Val-Borro et a l.|2010[ l, the baseline ripple on the Mars 
observations is rather large because of its strong continuum 
emission. This phenomenon was frequently observed by ground- 
based observations of planets. While in the cometary case, the 
baseline ripple has been removed with a polynomial fit, for Mars 
we determined the baseline frequencies by a normalized peri- 
odogram according to Lomb (1976) and subtracted 3 periods 
from the original spectrum. The influence of the line shapes on 
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Fig. 1. C ls O (left) and 13 CO (right) WBS spectra before baseline 
ripple removal. Both lines are in the same sideband. 



the determination of the baseline periods is negligible when tak- 
ing the large number of periods in the spectrum into account. The 
observed spectral lines have been modelled using a standard ra- 
diative transfer code: Mars was assumed as a perfect sphere sur- 
rounded by a set of hundred concentric atmospheric layers of 1 
km thickness each (as in |Rengel et al.|2008"] l. Within each layer 
the atmospheric temperature, pressure, and volume-mixing ratio 
of carbon monoxide have been assumed constant. The surface 
continuum emission was modelled as blackbody emission using 
a temperature distribution falling off towards the edge of the ap- 
parent disk according to T(a) = Tq X (1 - 0.2 x (1 - cos(a))), 
with a running from 0° (nadir) to 90° (limb) across the apparent 
disk (limb darkening, compare also Cavalie et al. 2008). 

The emission was obtained by integrating over the apparent 
disk using 64 concentric rings distributed unevenly over the disk 
and the limb region. The variation in the path lengths through 
the atmosphere were taken into account fully when calculating 
the radiation transfer for each ring. In our model the total contin- 
uum flux emitted by the surface depends purely on the choice of 
the temperature Tq, which defines the temperature scale for the 
temperature profile to be retrieved. We have adjusted To in such 
a way as to exactly match the total flux predicted by the 'Mars 
continuum model' of ~ 4230 Jy provided by Lellouch & Amri| 



(2008) and to match the temperature fall off towards the limb 
therein by a factor 0.2. The error of the modelled flux is 5%. 

Absorption coefficients for the CO spectral lines were calcu- 
lated using the HITRAN 2008 spectral line catalogue, keeping 
the terrestrial isotopi c ratios in it; i.e., C 16 Q/C 18 = 498.70 and 
12 CO/ 13 CO = 89.01. iKrasnopolsky et al] ( |2007| conclude that 
the deviations of these isotopic ratios on Mars compared to Earth 
are less than 2%. However, to account for carbon dioxide instead 
of air as broadening gas, the broadening parameters provided by 
the catalogue were multiplied by a factor of 1.4 according to 
|Nakazawa & Tanaka| ( fl982] i. 

To obtain a temperature profile and the mixing ratio of CO 
we employed Rodger's optimal estimation method (Rodgers 
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Fig. 2. The I3 CO and C 18 lines after removal of the baseline ripples. A radiative transfer model was fitted simultaneously to the 
spectra to retrieve temperature profile and volume mixing ratio of CO. The lower panels show the differences between model and 
observations. 



|1976|[1990| . This method uses so-called 'a priori' information - 
in our case the best estimate of the temperature profile and CO 
mixing ratio profile - which is then updated to allow for a best 
fit of the spectral line shapes (minimization of x 1 )- This a priori 
profile was taken from our MAOAM general circulation model 
(Hartogh et al. 2005) considering the exact observation date, ge- 
ometry, and time of the HIFI observations. The surface pressure 
averaged over the visible disk was 6.7 hPa. Compared to temper- 
ature profiles derived from the EMCD, we find temperature dif- 
ferences less than 3 K between 3-60 km. EMCD provides about 
5 K higher temperatures below 3 km and 3 to 5 K higher tem- 
peratures between 60-80 km. Up to 100 km, both models slowly 
merge to almost the same temperature. It is worth noting a high 
degree of coincidence of the averaged temperature profiles over 
the field of view of the telescope with the two models. Even 
though the one from EMCD represents monthly averaged fields, 
the one from MAOAM is based on an instantaneous snapshot, 
and generally the altitude-latitude distributions differ. This in- 
creases our degree of confidence in the simulated temperatures. 

A simultaneous fit of the two spectral lines allows retriev- 
ing temperature and mixing ratios independently, because of the 
rather different optical depth of the two lines, where the 13 CO 
line is optically thick (r = 6.3 in line centre) and the C 18 line 
being optically thin (t = 1.1). Figure [2] shows the simultane- 



ously fitted spectra of 13 CO and C ls O (and the residuals) after 
the removal of the baseline ripple. The retrieved CO mixing ratio 
amounts to 980± 100 ppm which is in agreement with the one de- 
tected by SPIRE observations during L s = 5° from 6 November 
2009 at 20:20 UT ( |Swinyard et~aT1[20T0| > under rather similar 
surface pressure conditions. The SPIRE value of 900 ppm has 
been used as an a priori input to the retrieval algorithm. Figure [3] 
shows the corresponding temperature profile and the averaging 
kernels. The latter provide information about the sensitivity of 
the retrieval versus altitude. Although the contribution of the a 
priori profile to the retrieved temperature profile is less than 10% 
below 60 km, it fits quite well to the profiles predicted by the 
GCMs. However, the differences are least with EMCD near the 
ground (~ 5 K) and with MAOAM near 60-70 km. Nevertheless, 
the observations provide about 12-15 K lower temperatures near 
65 km. The temperature inversion between 40-60 km predicted 
by the GCMs should be manifested in an emission feature in 
the centre of the CO lines, which obviously is not the case. The 
5% error of the model continuum flux translates into a roughly 
5% shift of the temperature profiles, i.e., all temperatures will be 
about 10 K higher or lower than the retrieved value. In the first 
case, the agreement with the model profile is better between 40- 
80 km (although still without temperature inversion) but worse 
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Fig. 3. Retrieved vertical profile of temperature (left) and averaging kernels (right). The averaging kernels express the response 
of the fitted temperature profile to a delta function perturbation in the true profile at a certain altitude (see Rodgers 1990 for a 
sound introduction of this quantity). Different colours are used to distinguish the averaging kernels belonging to perturbations at the 
different altitude levels. The blue temperature profile is derived from EMCD. 



below 40 km. In the second case the agreement is worse for all 
altitudes. 



4. Summary 

This work presents the first simultaneous retrievals of temper- 
ature and carbon monoxide in the Martian atmosphere derived 
from HIFI data. The temperature profile can be used as an in- 
put parameter for determining of concentrations of other gases 
observed by HIFI during the same period. Future work will in- 
clude all observed CO transitions in order to better constrain the 
temperature profile above 60 km and take advantage of the much 
wider opacity range for retrieving the vertical profile of CO. 
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